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ABSTRACT 



We use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as 
part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass 
of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and 
compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared 
to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). 
This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, 
for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, 
similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower 
(by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor 
gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from 
the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the 
possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z. 

Key words. Galaxies: evolution - Galaxies: high-redshift - Galaxies: ISM - Infrared: galaxies - Submillimeter: galaxies 
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1. Introduction 

Understanding the evolution of the dust properties and the dust 
content of galaxies through the cosmic epochs is crucial for con- 
straining galaxy evolutionary scenarios. The amount of dust has 
important implications for the star formation (SF). Dust allows 
the formation of low-mass stars in low-metalicity environments 
while it in hibits the formation o f massive stars, hence affects the 
IMF (e.g. lOmukai et al]l2005h . Dust greatly enhances the for- 
mation of many molecules whose transitions provide the main 
cooling mechanism for molecular clouds to form stars. Hence 
dust is one of the pre-requisites for enhanced SF. Dust also af- 
fects the detectability of high-z galaxies: dust extinction reduces 
the detectability of high-z galaxies in the rest-frame UV-optical 
bands, while dust thermal emission allows the detection of high- 
z galaxies at mm/far-IR wavelengths. The dust content can also 
be used as a proxy of the metalicity of the ISM, because refrac- 
tory elements in the ISM are mostly depleted into dust grains. 

The ESA Herschel Space Observatory (Pilbrattet al. 2010) 
offers the possibility to investigate the dust emission in high- 
z star-forming galaxies and in particular to measure their dust 
content by modeling the far-IR to submm spectral energy dis- 
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tribution (SED). We here present a first investigation of dust 
mass evolution by focusing on a sample of high-z submillime- 
ter galaxies (SMGs). They are very luminous (~ IO'-^Lq) (e.g., 
Hughes e t al. 1998 ; Pope et al. 2006), high redshift (z ~ 1-3.5) 
(Chapm an et al.l 120 05: Pope et al. 20 0^, gas rich and co mpact, 
massive galaxies ('Grev e et al.ii2005i: iTacconi et al.ll2008l) . Their 
star formation rates (SFR) are exceptionally high (~ 10^ Moyr"' ) 
and a re thought to c ontribute significantly to the cosmic SF at z ~ 
2-3 dChapm an et al. 2005). Thanks to their high luminosities 
they provide an excellent laboratory to investigate the evolution 
of dust in high-z galax ies. To mea sure their dust content we use 
data taken with PACS dPogUtsch e t al. 2010) as part of the PACS 
Evolutionary Probe guaranteed time key program combined with 
ground-based mm and submm data. By comparing the infeiTed 
dust masses with other properties of the SMGs, and with local 
galaxies, we infer important information on the physics of these 
high-z systems. We adopt the A-CDM concordance cosmologi- 
cal model (Ho = 70 km/s/Mpc, Qm = 0.3 and Qa = 0.7). 



2. Sample and data analysis 

We used deep 100 and 160 jum PACS GTO observations of 
the field GOODS-N and of the lensing cluster Abell2218. 
Observations, data reduction and source extraction are described 
in detail in Bertaetal. (2010). As mentioned above, we used 
a subset of PACS sources with submm detections taken from 
the literature. The SMG sample and counterparts association 
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are described in iMagnelli et alJ (l2010t) . The sample is restricted 
to sources with secure redshifts (all spectroscopic, except for 
three sources with good photo-z) derived from robust multi- 
wavelength associations. A few SMGs in the GOODS-N field 
have two optical, and eventually PACS, counterparts at the same 
redshift. For these we considered th e total system o f the two 
galaxies, which are likely interacting dPope et al .'2006'). by sum- 
ming their IR fluxes and stellar masses. We excluded three galax- 
ies (GN19, GN19b, GN30) with only two photometric points at 
A >100/im because their SED could not be properly constrained. 
Finally, we excluded a galaxy (azt23) with a peculiar SED, sus- 
pected to be dominated by a powerful AGN even at far-IR wave- 
lengths. We ended up with 12 SMGs in GOODS-N and 5 in 
Abell2218, three of the latter are different images of the same 
lensed galaxy. All our galaxies have PACS (100 and 160;um) 
and SCUBA 850yum photometry, and nine also have SCUBA 
450yum and/or AzTEC 1.1mm and/or MAMBO 1.2mm detec- 
tions. Redshifts range from 0.5 to 4, with a median value of 2. 

Both to expand the sample and to check that our dust 
masses are not subject to systematics associated with the spe- 
ci fic set of dat a, we also included six SMG galaxies observed 
by iKovacs etalJ (12006) at 350//m with SharcII, which have at 
least two additional (sub)mm photometric poi nts, and one MIPS 
(both 70 and 160jum) detected SMG from the lYan et all (l2010l) 
sample, with also a MAMBO 1 .2mm observation. 

Dust masses were derived by fitting and normalizing the 
100/im-to-(su b)mm photometr y to the SED library generated 
by GRASIL dSilvaet al.lll998h . This is a chemospectrophoto- 
metric code that uses realistic and physically grounded stel- 
lar and dust distributions and geometries and performs the full 
radiative transfer calculations, eventually providing the com- 
plete SED of various galaxy models. Recently, GRASIL has 
been convolved w i th models of dust evolution in galaxies 
dCalura et all 120081: ISchurer et all l2009h . The model carefully 
takes into account the different emissivities of different popu- 
lations of dust grains, which also depend on their size, com- 
position, and temperature. This effectively results in an average 
emissivity index B ~ 2, cons istent with va rious prev ious studies 
dSilva et al.lll998llDraine &" Li 2007; Cle ments et al . 2010). The 



dust masses in ferred by us broadly agree with those obtained by 
iMichalowski et al.. (2010a), who also use GRASIL, although the 
lack of data at 24/im < A < 850//m in most objects prevents 
them from tightly constraining the far-IR bump. 

The inferred dust mass is compared with the stellar mass, gas 
mass, and gas metalicity as inferred for a subset of the objects in 
our total sample by ancillary data presented in the Appendix. We 
also compared the dust properties of SMGs with those of local 
star-forming galaxies, both spirals and ULIRGs, whose proper- 
ties are reported in the Appendix. Here we only mention that 
the methods to derive the dust mass, stellar mass, gas mass, and 
metalicity are the same for all samples, so that their comparison 
is done in a fully consistent way. 

3. Results and discussion 

Figure [T] shows the stellar mass vs the dust mass for high-z 
SMGs (black circles), local ULIRGs (blue squares) and local 
spirals (red triangles). The dust masses estimated for local spirals 
nicel y agree with those ob tained through the very detailed analy- 
sis of iDraine et al.l (l2007l) of the same sample. Most of them also 
agree fairly well with the evolutionary pattern for spiral galax- 
ies predicted by the model of Caluraetal. (2008) (cyan solid 
line). The dust masses were compared with those obtained by the 
more common assumption of a single temperature dust distribu- 
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Fig. 1. Stellar mass versus dust mass. Blue squares show local 
ULIRGs, red triangles refer to local spirals, while black circles 
correspond to high-z SMGs. Large gray, magenta and green cir- 
cles mark respectively photo-z, the triple image in Abell2218 
and blended SMG systems. Median Icr error bars for the diff'er- 
ent samples (same color code) are shown on the left. The solid 
cyan and dashed (dot-d ashed) orange lines show the predictions 
of Calur a et aP (l2008h model for spirals and proto-ellipticals 
with mass of 10" {\0^^)Mq. 



tion. Assuming an emissivity index /3 - 2 (for consistency with 
GRASIL) and a single, absolute average emissivit y of dust at 
125/[/mof2.64(DunneeLan^2003")orl.87m-kg-i (ILi & Draind 
2001), dust masses are on average 0.3 or 0.15 dex, respectively, 
lower than those in ferred through GRASI L (lOOyum flux not used 
in the fitting as in IMagnelli et aLluOl O*). Independently of the 
adopted method, there are clear, systematic diff'erences in terms 
of dust content for the different samples. The orange dashed and 
dot-dashed lines in Fig. [1] sho w the expected evolution of proto- 
ellipticals (ICalura et al.l 120081) with different mass. These mod- 
els predict a larger dust content during the active phase of ellip- 
tical formation compared to normal spirals. Nonetheless, they 
can hardly reproduce the large Mjust/Mstar observed in most 
SMGs. The dust-to-star excess in SMGs emerges more clearly 
in Fig. |2^, which shows the distribution of Mdust/Mstai for the 
different samples. Local ULIRGs have a larger dust content than 
spirals for a given stellar mass, with Mdust/Mstai - a factor of ~4 
higher on average. This result is in contrast with IClements et aTl 
(120101) . who claimed that the dust content in ULIRGs can be sim- 
ply explained with the combination of two spirals (the different 
result probably arises because Clements et al. use only K-band 
photometry as a proxy of the stellar mass, while we perform 
the full SED fitting). Submillimeter galaxies are much more ex- 
treme: their Mdust/Mstai- is on average higher by a factor of ~30 
than in local spirals, and by a factor of ~6 than in local ULIRGs. 
This large dust content is difficult to account for with the galaxy 
evolutionary model of Calura et al] d2008l) . 

High-z SMGs are known to be very gas ric h, with gas 
fractions approaching 50% dTacconi et al.l l2008h . The high 
Mdust/Mstai- values suggest that the large gas masses hosted 
by SMGs have already been highly enriched with metals and 
dust; i.e. the excess of dust mass in SMGs is not a conse- 
quence of SF episodes following recent inflows of metal-poor 
gas accreted during mergers (in contrast with the expectation 
of some models, Montuori et al. 2010). This can be tested 
by measuring the dust-to-gas ratio for those SMGs with CO 
millimetric data. By assuming the same CO-to-Hi conversion 
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Fig. 2. Panels a), b), d): Mdust/M.tai-, Mdust/Mn, and Mdust/Mz.gas 
distributions, respectively, for high-z SMGs (black hatched), lo- 
cal spirals (SINGS, red open) and local ULIRGs (blue shaded). 
Mhj is computed adopting a = 4.3. Adopting » = 1 moves 
the black and blue histograms by the amount indicated by the ar- 
rows. Panel c): metal abundances for the sample of high-z SMGs 
of Swinbank et al. (2004) (black hatched histogram) and that of 
local ULIRGs of Rupke et al.. (.2008.) (blue shaded histogram). 
Submillimeter galaxies in common with our sample are marked 
by arrows. The red solid and green dashed lines show the met- 
alicity distributio n for local spirals (S PSS) and z ~ 2.2 star- 
forming galaxies (iMannucci et al. 2010h. respectively, with stel- 
lar masses similar to those of our SMGs. 



factor a - 4.3 M0(K km s"' pc^)"' for all samples (typ- 
ically applied to normal spirals) the resulting distribution of 
Mdust/Mnj is shown in Fig. |2j), similar for the three samples. 
If one adopts for ULIRGs and SMGs the conversion factor 
a ~ I M0(K km s"' pc^)"' (thought to be more appropri- 
ate for these classes of objects, ISolomon & Vanden Boutll2005t 
iTacconiet all 120081) . then ULIRGs and SMGs have Mdust/Mn, 
significantly larger than local spirals, as indicated by the arrows. 
It is not possible to infer the total gas mass for SMGs by in- 
cluding the fraction of HI, because this is not observed; but sim- 
ilarly to local UL IRGs, the bulk of the ga s is probably in the 
molecular phase (ISanders & Mirabel! 1 1996 1) in contrast to spi- 
rals, which may have a considerable fraction of gas in the atomic 
phase. As a consequence, by including Mhi to obtain the dust-to- 
total gas ratio, the resulting Mdust/Mgas distribution of ULIRGs 
and SMGs would be further skewed towards higher values rel- 
ative to spirals. A dust-to-gas ratio higher than in spirals has 

region of loc 
^I2008h . It 

sibly be ascribed to enhanced dust production by SNe and by 
the first generation of AG B stars in young stellar systems (e.g. 
iMichalowski et al]l2010 b'). and/or to higher accretion of metals 
on dust grains in these dense environments. In either cases, the 
finding that the dust-to-gas ratio in SMGs is not lower than in 
local spirals supports the interpretation of their "dust richness" 
(relative to the stellar mass) in terms of "gas richness", possi- 
bly along with a higher Mdust/Me,- Summarizing, the high dust 
content observed in SMGs mirrors their high gas content, the lat- 
ter with a dust-to-gas ratio similar or higher than normal spiral 
galaxies. 



also been observed in the centr al (dense) region of local lumi- 
nous infrared galaxies (LIRGs) ( Wilson et al.ll2008l) . It can pos- 



The latter result is in contrast with the low gas metalicity 
observed in SMGs though. Indeed, based on their optical neb- 
ular lines SMGs have metalicities similar or lower than other 
strongly s tar- forming system s, such as ULIRGs in the local 
Universe (iRupke et al. 20081). as well as Lyman break galaxies 
(LB Gs), BzKs (seelDaddi et al.ll2004l) and ULIRGs at high red - 



shift (lErb et al.l l2006' 'Mannucci et al.ll2O10HCaputi et alEOOSh . 
which in turn are lower than local, moderately star-forming spi- 
rals with the same stellar mass. This is illustrated in Fig. |2};, 
where the metalicity distribution of different classes of ob- 
jects was measured with the same method (based on optical 
nebular lines) and adopting the same metalicity scale as dis- 
cussed in the Appendix, and within similar stellar mass inter- 
val (M* ~ 10"Mo). For the SMGs we explo ited the sample 
with rest-frame optical spectroscopy given in ISwinbank et al.l 
(2004), which unfortunately overlaps with our sample only for 
three objects (marked with black arrows). That galaxies with 
enhanced star formation are characte rized by lower metalic- 
ity has been noted by various author s (iKewlev & Ellison 20081: 
'Rupke et al.l l2008l: ICaputi et al.ll2008h . Recentlv IMannucci et al.l 
(2010) showed that the anticorrelation between metalicity and 
SFR actually accounts for most of the evolution of the mass- 
metalicity relation at high redshift: high-z galaxies (at least up 
to z~2.5) appear more metal-poor simply because characterized 
by higher SFRs as a consequence of selection effects. This anti- 
correlation between SFR and metalicity is ascribed to inflow of 
metal-poor gas that boosts SF in galaxies. However, the inferred 
"low" gas metalicities in SMGs are in striking contrast with the 
high dust content inferred from their far-IR emission. 

In Fig. [3] we further investigate this discrepancy by compar- 
ing the dust-to-gas ratio with the gas metalicity. Red triangles are 
all local galaxies for which the complete required information is 
available; most of the upper limits are due to the lack of detec- 
tion or of data for the molecular gas. The dotted line shows the 
trend expected by assuming that the dust content scales linearly 
with metalicity following the relation (Draine et al. 2007) 

Mdust/Mgas * 0.01 X (0/H)/(0/H)mw. (1) 

Actuallv. lDraine et al.l(l2007h and lHunt et"an (l2005h have shown 
(on a wider sample) that the relation traced by real galaxies is 
much steeper than expected by Eq.[T] as also hinted by the (much 
smaller) subset of local spirals/dwarfs shown in Fig. |3] Local 
ULIRGs (blue symbols) tend to have higher dust-to-gas ratios 
for a given metalicity, the deviation depending on whether the 
assumed CO-to-H2 conversion factor is or = I (open squares) or 
a - 4.3 (crosses). In the case of SMGs (black symbols) there are 
only two objects for which the complete required information 
is available (whose gas metalicity is actually an upper limit), 
shown both in the case of c = 1 (solid squares) and a - 4.3 
(solid circles). However, it is clear that their dust content is far 
higher than expected from their metalicity. The HI mass is not 
included for ULIRGs and SMG s, but the effect is very small 
(0. 1 dex, ISanders & Mi rabel 1996). The orange and cyan tracks 
in Fig. [3] show the expec ted evolution of the dust content and 
metalicity as predicted bv lCalura et al.l (l2008h for ellipticals and 
spirals, where it is evident that the large discrepancy observed 
in SMGs cannot be accounted for by models. In short, in SMGs 
the mass of metals inferred from the dust content is much larger 
than the mass of metals inferred from the gas phase. 

To further quantify this discrepancy, in Fig.|2}J we show the 
distribution of Mdust/Mz.gas, where Mz,gas is the mass of metals 
in the gas phase (i.e. Z x Mgas). We assumed for all objects a 
conversion factor a - 4.3, the arrows showing the effect of using 
Of = 1 for ULIRGs and SMGs. Clearly in SMGs Mdust/Mz,gas 
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is over one order of magnitude higher than in spirals and well 
above unity. Local ULIRGs share a similar problem, but not as 
extreme as SMGs. When computing Mz.gas, we assume the same 
metal abundance for the neutral and ionized gas phases. This 
can be false, as the former may be ch aracterized by lower (by a 
factor of ~7, iLebouteiller et al.l2009h abundances than the latter. 
However, this effect does not change our results. Indeed, dwarf 
low-metalicity galaxies (see the Appendix), where most of the 
gas is in the atomic phase, would move towards higher values of 
Mdust/Mz.gas by 0.85 dex, while normal spiral galaxies, whose 
atomic gas fraction is around 50%, would only move by 0.24 
dex. The ratio Mjust/Mz.gas in spiral galaxies and dwarfs would 
still be lower than in SMGs, where the bulk of the gas is expected 
to be in the molecular phase (as discussed above). 

There are some possible scenarios that could explain this 
puzzling result. One possibility is that dust masses are erro- 
neously estimated because the dust properties in SMGs are very 
different from those assumed in our models. The average dust 
emissivity may be higher. However, this difference would have 
to be quite dramatic, because the dust masses are off by more 
than an order of magnitude with respect to those expected by the 
metaliciti es. In ULIRGs (extreme starbursts) the emissivity in- 
ferred bv [elements et al.l (12010.) is similar to that assumed here. 

Another possibility is that as mentioned above the high den- 
sity in SMGs favors the growth of dust out of metals in the 
ISM, which also decreases the effective metalicity in the gas 
phase. However, the standard dust depletion factors in the dif- 
fuse ISM are already high, with most metals already locked into 
dust. Hence, for a given metalicity the dust mass cannot grow by 
more than a factor of about two through this mechanism. 

Finally, the disagreement between the dust content and the 
gas metalicity can be caused by dust obscuration. Likely, the 
bulk of the gas in SMGs is metal rich. The associated dust rich- 
ness and the compact configuration make most of the ISM op- 
tically thick at visual wavelengths. As a consequence, the op- 
tical nebular lines used to infer the gas metalicity only probe 
the outer parts of the star-forming regions, which are probably 
more metal-poor. Paradoxically, the higher metalicity of these 
dense systems and their associated high dust content may pro- 
duce the apparent effect of metal-poorness when observed at op- 
tical wavelengths. T his scenario is supported by the findings of 
ISantini et alJ (l2009h . who showed that luminous IR galaxies have 
higher mid-IR-based SFR estimates compared to those obtained 
by correcting the optical-UV light. 

If confirmed, the latter two scenarios would not support the 
current interpretation of the low metalicity observed in ULIRGs 
and SMGs in t erms of inflow of meta l-poor gas boosting the star 
formation (e.g. lMontuori et al]l2010l) . More generally, if the dis- 
crepancy between gas and dust metal mass also applies to other 
populations of high-z starburst galaxies (e.g. BzKs and LBGs) 
this would prompt a general revision of the interpretation of the 
apparently low metalicities observed in high-z systems. 
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Appendix B: Ancillary data 

Here we first describe the ancillary data used for SMGs (stellar 
mass, gas mass and metalicity) and afterwards the local samples 
used for comparison. 

B. 1 . Submillimeter galaxies 

Stellar masses of SMGs are computed by fitting the optical- 
to-near-IR photometry compiled b y the PEP Team to synthetic 
models, which is fully described in lFontanaet al.l (|2006), fixing 
the redshifts to tho se of the SMG s. For consistency with the the- 
oretical models of ICalura et al.l (2008) with which we compare 
our results, we adopt a Scalo IMF, which implies stellar masses 
higher than those calculated assuming a Salpeter IMF by a factor 
1.15. 

When possible, we obtained information on the gas mass. 
For the SMG sample we were able to find CO millimetric data 
for 12 galaxies (the three measurements of the tri ple image of the 
lensed galaxy in Abell2218 are counted once) (iBothwell et al.l 
20091: iKneib et al.1 120051: iDaddi et al.1 120091: iGreve et al.1 120051: 



Knudsen et al.l 1200 9). To derive the molecular gas masses we 



followed the prescriptions given in the associated references, but 
we normalized them to the same CO-to-H2 conversion factor. In 
particular we show the results assuming (conservatively) a CO- 
to-H2 conversion factor a = 4.3 M0(K km s"' pc^) ', which 
applies to normal spirals; but we also show the effect of as- 
suming ff = 1 Mq(K km s"' pc^)"', w hich is thought to be 
more appropriate fo r ULIRGs and SMGs dSolomon et al.ll 19971: 
lTacconietall2008h . 



We inferred the gas metalicities through the optical nebular 
line emission by using the strong-line methods. Within this con- 
text, and especially when comparing different samples, it is most 
important to use consistent (inter-)calibrations and, if possible, 
use the same diagnostic em ission line ra tios. Here we adopted 
the calibration s obtained by iNagao et al.l (2006) and refined by 
Ma iolino et al.l 02008). who inter-calibrate all strong line diag- 
nostics to the same scale. For the SM Gs we used as m e talicity 
diagnostic the ratio [NII]/Ha given in ISwinbank et al.l ( |2004|) . 
An AGN may easily aff'ect the nebular line ratios and make 
them deviate from the metalicity calibration for star-forming 
galaxies. For this reason we removed targets classified as AGNs 
in S winbank et al. (2004). We also removed one object with 
[NIIl/Ho- =0.65, because such a high ratio is only observed in 
AGNs or shocked regions. In total, we obtain gas metalicities for 
a sample of 10 SMGs, three of which also belong to our SMG 
sample with dust masses. The [NIIl/Ha ratio of these objects 
may still be aff'ected by an unrecognized AGN, but this would 
make our results conservative, because the contribution by an 
unrecognized AGN would make the gas metalicity apparently 
higher than its actual value. 

B.2. Local galaxies 

To compare the dust content in high redshift SMGs with nor- 
mal and actively star-forming galaxies in the local Universe, we 
assembled two sets of data: 

- a sample of 26 loca l spirals out of t he SINGS sample 
dKennicutt et al.1 120031; iDale et al.ll2007t iDraine e"tani2007l) 
with full multi-wavelength photometry and also with submm 
data for a proper dust mass measurement; 

- a sample of 24 local ULIRGs taken from the lClements et al.) 
(2010) sample with at least three submm photometric points. 
Optical and near-IR fluxes needed to compute the stellar 
mass of these galaxies where collected from public archives. 

Stellar masses were inferred by the same method and IMF as 
SMGs. 

For the l ocal spirals, we took the HI and H2 mass es- 
timates from iKennicutt elaP (l2003l) . w hile for the ULIRG s 
we collected H2 rnass estimates from Sanders et al.l (Il988h . 
'Sanders et al." (1991"), 'Solomon et dj (ll997h . lEvans et al.l(l2002l ) 
and PacadoDoulos et al. (2010) . Here we also assumed a CO- 
to-H2 conversion factor a = 4.3M0(K km s"' pc^)"', but 
for ULIRGs we also show the eff'ect of assuming a = 
lMa(Kkms-ipc2)-i. 

Concerning the gas metalicity in the local spiral sample, 
there are only a few objects (six after removing possible AGN- 
dominated galaxies) out of the SINGS sample with submm data 
that also have integrated spectroscopy dMoustakas & KennicuttI 
2007; Kennicutt 1992), which could be used to infer the gas met- 
alicity with the same strong line method and calibration scale 
adopted for SMGs. For these objects we also exploited the di- 
agnostics involving [OII]3727, [OIII]5007 and H/3, which help 
to better constrain the metalicity along with [NII]/Hq', while 
they are still fully consistent with th e global inter - calib rations 
presented in Mai olino et alJ (l2008h . iDraine et all (12007:) uses 
the metalicities for a more extended set of SINGS galaxies, 
but by adopting different calibrations, which therefore cannot 
be used for direct comparison with the SMGs. To expand the 
sample of objects with dust masses and metalicity (especially 
in the low-metalicity range), we included nine dwarf and spi- 
ral galaxies of the SINGS sample without submm photometry 
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(iDraine et a 



iDraine etal 



2007 ). There the Mdust estimates were taken from 
(I2007h (whose dust masses are consistent with ours 
for the objects in common, as discussed in the text), and their 
integrated emission line fluxes, which were required to place 
them in the same metalicity scale as SMGs a n d ULIRGs, were 
collect ed from Moustakas & KennicuttI (l2007h . iMTUer & Hodg3 
(11996") and'Miller (,1996l) . 

For the ULIR Gs we collected spectroscopic data from 
iRupke et alj (l2008h . We removed galaxies classified as AGNs 
based on the diagnostic diagrams. 

In order to have a more comprehensive statistic (only 6 out 
of the 26 local spirals from the SINGS sample have integrated 
spectroscopy which can be used to compute their metal abun- 
dance), the metalicity distribution of local spirals shown by the 
red, solid li ne in Fi g. [2}: is obtained from the large SDSS sample 
presented in lKewlev & Ellison (2008) in the same mass range as 
SMGs and ULIRGs, with th e same metalicity calibrations and 
scale used in lMaiolino et al.l (2008), i.e. consistent with the other 
samples. The metalicity distribution of LBGs and BzKs at z~2. 2 
(with logM* ~ IO.6M0) is taken from 'Mannucci et al.' ('2010^ 
who adopt the same scale and calibrations as .MaioUno et al] 
(l2008h . 



